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' Abstract. Using the method of Darboux transformations (or equivalently supersymmet- 

' ric quantum mechanics) we obtain an explicit expression for the propagator for the one- 

dimensional Sehrodinger equation with a multi-soliton potential. 



Key words: Darboux transformations; SUSY QM; soliton potentials; propagator 
2000 Mathematics Subject Classification: 81Q05; 81Q60 



^ , 1 Introduction 

> 

00 . 

. It is nowadays incontestable that ideas of supersymmetry play an important role in modern 

CN I theoretical and mathematical physics. Super symmetric quantum mechanics (SUSY QM), intro- 

duced by Witten ^ as a toy model in supersymmetric quantum field theory, is a very efficient 
l/^ . tool for studying different properties of non-relativistic quantum systems. In particular, as it 

I was recently shown |2] it may become an essential ingredient of complex quantum mechanics 

which is currently under development since it may "cure" such a "disease" of non-Hermitian 
Hamiltonians as their non-diagonalizability, and can remove spectral singularities from the conti- 
I nuous part of the spectrum. On the other hand, it is well known j^l that SUSY QM is basically 

' equivalent to the method of Darboux transformations [1] well known in connection with soliton 

^ theory [H]. We mean by soliton potentials such potentials of the one-dimensional stationary 

Sehrodinger equation from which soliton solutions of the Korteweg-de Vries equation may be 
constructed [H] and which are SUSY partners of the zero potential. They find an application in 
^ ' quantum field theory for describing processes where solitons may play an essential role [HI. 

It is not an exaggeration to say that everything in nature undergoes time-evolution. There- 
fore, one of the crucial questions of any physical theory is to describe how a physical phenomenon 
evolves with time. In quantum mechanics time dependence of the wave function may be de- 
scribed with the help of the propagator which is nothing but the evolution operator in coordinate 
representation. As far as we know, the first attempt to find the propagator for a multi-soliton 
potential is due to Jauslin [Jj. His method is based on an integral formula which relates so- 
lutions of two Sehrodinger equations whose Hamiltonians are SUSY partners. Unfortunately, 
this approach when applied to the Sehrodinger equation leads to divergent integrals. There- 
fore the author found the propagator for the heat equation with the one-soliton potential. The 
Sehrodinger equation may be considered as the heat equation with the imaginary time. In this 
respect the following question arises: whether or not the Jauslin's result after the replacement 
t ^ it gives the propagator for the one-soliton potential? We want to stress that the answer to 
this question is not trivial since such a replacement at the level of the Jauslin's integral trans- 
formation leads to divergent integrals. Our analysis shows that the answer to this question is 
positive. 

We would like to note that the problem of finding propagators is more difficult than solving 
the stationary Sehrodinger equation for a time-independent Hamiltonian, since knowledge of 
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the propagator allows us to solve the Cauchy problem for the non-stationary equation with the 
same Hamiltonian, but with an arbitrary initial condition. This is similar to the Green function 
of the stationary equation; knowing it allows us to write down a solution of the inhomogeneous 
equation with an arbitrary inhomogeneity. Recently a method has been proposed for finding the 
Green function for a SUSY partner Hamiltonian [5]. In this paper, using a particular example 
of soliton potentials, we show that the method of SUSY QM is very helpful for finding exact 
propagators for Hamiltonians related by SUSY (or equivalently Darboux) transformations. We 
would like to stress that although many aspects of multi-soliton potentials are well-studied in 
the literature a closed form of the propagator as far as we know has not yet been published. 

To make the paper self-contained, first we briefiy review the method we are using and give 
a general expression for the propagator (Section 2 and the beginning of Section 3), and then 
apply it for finding the propagator for the multi-soliton potential (Section 3). Finally, we apply 
our general result to the case of one- and two-soliton potentials (Section 4) and suggest a formula 
for A^-soliton potential. 



2 Darboux transformations and soliton potentials 

In this section we review briefiy the method of Darboux (SUSY) transformations El Also 
we construct a differential (i.e. Darboux) transformation operator for obtaining solutions of the 
Schrodinger equation with a soliton potential from solutions of the free particle equation we 
need in the following sections. 

Consider two one-dimensional Schrodinger equations 

(/lo - ^)V^s(x) = 0, ho = -d'', + Vo{x), (1) 
{hi - E)ipEix) = 0, hi = -dl + Viix). (2) 

Suppose one knows solutions to equation (^. Then solutions to equation Q can be found 
by acting with the transformation operator (we denote it as L) on solutions to equation 
(fiE = LipE- The main relation defining L is the intertwining relation 

L/io = hiL. (3) 

An essential point of the method is the choice of the operator L. If we restrict L to be 
a differential operator it becomes Darboux transformation operator (see e.g. (HI). In this case 
the potential Vi cannot be arbitrary, and should be found together with the transformation 
operator from the intertwining relation 12)). If L is a first order differential operator, the result 
is well-known (see e.g. pi^ITU]): 

L = -d^ + w{x), Vi = Vo- 2w'{x), 

where 

w(x) = (Inn)a;, hou = au. 

Function u = u{x) and parameter a are called transformation function and factorization constant 
respectively. 

Since the procedure is independent on the initial Hamiltonian, it can be repeated as many 
times as one desires. So, we can take hi as the initial Hamiltonian for the next transformation 
step to get /i2 and so on. In this way one arrives at chains of transformations. It is remarkable 
that the resulting action of a chain may be expressed in terms of solutions of the initial equation 
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only. If all factorization constants are different from each other, one obtains Crum-Krein jll) 
formulas 

ipE = LtpE = W{ui,U2, . . . ,un,iPe)/W{ui,U2, ■ ■ ■,Un), (4) 
Vn = Vo-2 [log W{ui,U2, un)]" ■ 

Here and in the following the symbol W denotes a Wronskian, and uj = Uj{x) is an eigenfunction 
of the initial Hamiltonian (transformation function) 

hQUj{x) = ajUj{x), j = 1,2, . . . , N. 

To get soliton potentials one starts with the zero initial potential, Vo(x) = 0, and uses the 
following set of N (which is supposed to be even) transformation functions ^ El E] 

U2j~i{x) = cosh(a2j-i2; + b2j^i), (5) 
U2j{x) = sinh(a2jX + b2j), j = 1,2,. . . , N/2. (6) 

They are solutions to the Schrodinger equation with the zero potential corresponding to eigen- 
values Ej = —Oj < 0, which are just the points of the discrete spectrum of = —d"^ + Vn{x). 
An orthonormal set of its discrete spectrum eigenfunctions is given by 

1/2 

^n[x)=\— a„-a. — ^ . 7) 

Here W^'^\ui,U2, ■ ■ ■ ,Un-i,Un+i, ■ ■ ■ ,un) is the Wronskian of order — 1 obtained from the 
Wronskian W{ui,U2, . . . ,uj\f) by dropping the function 

We also need the continuous spectrum eigenfunctions of /iat, which should be found by acting 
with the operator L ^ on plane waves ipk{x) = l/\/27r exp(— ifcx). A: E M 

V^k{x) = I Lipkix), (8) 

E = , ak = —clI, k = 1, . . . , N. 

The set of functions {(pnix), n = 1, . . . , N} and {ipk{x), k £ M} forms a complete and orthonor- 
mal set in the Hilbert space of square integrable functions on the whole real line. 

It is interesting to note that for particular values of the parameters aj a multi-soliton potential 
may have a shape of a multi-well potential thus presenting an example of a multi-well exactly 
solvable potential. 




3 Propagator for a multi-soliton potential 

We use the definition of the propagator K{x,y;t' ,t") of the Schrodinger equation as the co- 
ordinate representation of the evolution operator (see e.g. 23|)- If one knows solutions of the 
non-stationary Schrodinger equation with a time-independent potential, which form a complete 
and orthonormal set (for continuous subset, if present, normalization is understood in the sense 
of generalized functions) in the Hilbert space of square integrable functions on the whole real 
line, then the propagator is given by 

K{x,y;t',t")=9{t'-t") 



N 



^i^n(.X,t')'ll^n{y, 



t") + 



.n=0 



dkipi{x,t')My,i") 
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Here {tpnixjt)} is the discrete part of the basis and {ipkixji), k € M} is the continuous one. 
Everywhere we shall assume t' > t" and drop the step function 6{t' — t"). 
According to this formula, the propagator has two contributions 

K{x,y;t',t") = Kd{x,y-t',t") + K,{x,y;t',t"). (9) 

The first term at the right hand side of formula Q, K^, corresponds to the discrete part of the 
basis, and the second one, Kc, proceeds from the continuous part. Our crucial observation is that 
the method of Darboux transformations gives both the discrete and the continuous parts of the 
basis for the transformed equation if a complete set of eigenfunctions for the initial Hamiltonian 
is known. 

For soliton potentials the discrete part is given by and the continuous one is given by (jEJ. 
Therefore 

KdN = Y,^*n{x,t')ipn{v,t"), K,N= / dkipl{x ,t')ipk{v ,t") , (10) 

n=0 

where the additional subscript N labels the order of the Darboux transformation, which in our 
case coincides with the number of solitons. 

We compute first the value Kcn- After inserting (pk{x) from into (jlOj) we interchange 
derivatives with the integrals over the momentum k. This allows us to present the contribution 
from the continuous spectrum as an action of the transformation operator L on an integral 



KMX, y,t,t}- L^Ly dk ^ ^2)(^2 + „2) . . . (^2 + 

roD ^ik{x-y)+ik^{t'-t") 



L^Ly / dk—-r- — g — 2" — 2—. (11) 



27r " (fc2 + a2)(^2+„2)...(^2 + 



N) 



Here is given by @ and Ly is obtained from L,j. by the replacement x — > y. The integral 
in ifTTj) is calculated in Appendix. Thus, the contribution from the continuous part of the basis 
has the form: 

n=l \ m^n 

Here t = t' - t" and 

dp—o ^ = -e^" Re[e4^^'^erfc {aV2 + ixV2)\. (13) 

+ a 

To compute KdN, we simply replace in (fTUI) the discrete basis eigenfunctions ipi according 
to O and take into account their time dependence, which yields 

K fr ,2 2,^ W-(")(X)W-W(y) 

K<iN-^\^ 11 \^n-a,\\ ) e t-t-t. (14) 



Thus, we see that the propagator for the multi-soliton potential is expressed in terms of the 
error function, derivative of the error function and solutions corresponding to the discrete part 
of the basis. In the next Section, we suggest a simpler formula for the propagator. 
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4 Particular cases 

Let us consider how expressions ()12() and (|14|) lead to the propagators for the one- and two-soliton 
potentials. 

The one-soliton potential 

V{x) = -2aVcosh2(ax) 

may be obtained from the zero potential with help of the first order Darboux transformation 
with the transformation function u{x) = cosh(aj;). The only discrete spectrum eigenfunction is 
ipo{x) = (a/\/2)/ cosh(ax). 

After some simple calculations we obtain from ©, (fT ^ -(fT Hl the propagator for the Schro- 
dinger equation with the one-soliton potential 

Ki{x,y-t",t') = -^=e^-^ + -P-^ — ^ [erf + (a) +erf _(a)] , (15) 

VAttU 4cosh(ax)cosh(ay) 

t = t'- t". 



Here erf -|-(a) is the error function taken at a special value of the argument 
erf -|-(a) = erf 



aVii ± - y) 



Now two observations are in order, (i) We state that up to the replacement t ^ it our result 
is in the perfect agreement with that of Jauslin [Jj. This fact justifies the Jauslin's procedure of 
regularization of corresponding divergent integral. And (ii) we see that the propagator H15() has 
two contributions. The first term coincides with the propagator for the free particle and, hence, 
the second one describes just the one-soliton perturbation of the free particle at the level of the 
propagator. 

In case of the two-soliton potential the method of Jauslin becomes very involved, since it is 
not clear how corresponding integral may be calculated. In contrast, our formulas (|12|) and (|14)) 
give an explicit expression for the propagator. To use formula ()14() . we first calculate necessary 
Wronskians with the transformation functions where for simplicity we choose 6i = 62 = 0, 
i.e. ui{x) = cosh(aix), U2{x) = sinh(a2a;) and then after some algebra we get the propagator for 
the two-soliton potential 



K2ix,y;t",t') 



1 



le 4t 



01(02 — of) sinh(a2x) sinh(o2y)e*"i* 

AW{x)W{y) 
02(0^ — of) cosh(a23;) cosh(o2y)e*"2* 



+ 



+ 



AW{x)W{y) 



[erf _|.(oi) + erf _(oi)] 
[erf + (02) + erf _ (02)] , 



(16) 



t = t'- t". 



Once again we see that the propagator for the two-soliton potential has a structure similar to 
that of the one-soliton potential. Both the propagator ()15() and H16|) have two contributions: the 
exponential part which is the free particle propagator and the part responsible for the soliton 
perturbation of the zero potential. 

This result suggests us the form of the propagator of the general A^-soliton potential as a sum 
of two contributions, namely the free particle propagator in the form of the exponential function 
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and an A^-soliton perturbation: 



KN{x,y;t",t') 




VI7(")(3;)VF(")(2/) 



1 i{x-v)'^ 



W{x)W{y) ^ 





On 

4 



[erf + (a„) + erf _(afc)] 



n 



N 




t = t'- t". 



(17) 



We have checked this formula for = 3,4 by the direct substitution into the Schrodinger 
equation. 

5 Conclusion 

In this paper using the known approach based on supersymmetric quantum mechanics (or equi- 
valently, the method of Darboux transformation) we construct propagators for multi-sohton 
potentials. In particular, we present explicit expressions for the case of one- and two-soliton 
potentials. While comparing these results with the ones previously published by Jauslin [7], 
we conclude that the transformation t ^ it (Whick rotation) is justified for regularization of 
divergent integrals of a special type. Finally, we suggest an explicit expression for the propagator 
for the general multi-soliton potential. 
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A Appendix 

Here we calculate the integral 

/■oo ikx+ik'^t 




The identity 




reduces it to a sum of simpler integrals of the form: 




One can find the value of the integral at the right hand side of this equation in jl4j 




dp- = -i7re2"("-2«)erfc \V2ia - ix 

p + ia 



which leads just to formula (|13jl for the function 
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